Unitary quantum gates: a measure theoretic approach 
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Non-local properties of ensembles of quantum gates induced by the Haar measure on the unitary 
group are investigated. We analyze the entropy of entanglement of a unitary matrix U equal to the 
Shannon entropy of the vector of singular values of the reshuffled matrix. Averaging the entropy 
over the Haar measure on U(N 2 ) we find its asymptotic behaviour. For two-qubit quantum gates 
we derive the induced probability distribution of the interaction content and show that the relative 
volume of the set of perfect entanglers reads 8/3-7T ~ 0.85. We establish explicit conditions under 
which a given one-qubit bistochastic map is unistochastic, so it can be obtained by partial trace 
over a one-qubit environment initially prepared in the maximally mixed state. 
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I. INTRODUCTION 

Unitary quantum gates form key ingredients of any 
quantum algorithm, so they are widely used in the the- 
ory of quantum information A unitary gate acting 
on a bipartite system A <8>B is called local, if the unitary 
matrix has a form of the tensor product, U = V a ® Ub- 
To process quantum information between both subsys- 
tems one needs to use non-local gates, which arc not of 
the product form. 

Classification of unitary quantum gates is a subject of 
a considerable interest 2-5]. Quantification of the non- 
local properties of unitary gates has been initiated by 
Zanardi and co-workers [6|— lS| , while several other mea- 
sures of non-locality were introduced and analyzed in a 
seminal paper of Nielsen et al. Q. Since local unitary 
gates cannot produce quantum entanglement, the non- 
local properties of a given gate U may be characterized 
by the average (or maximal) degree of entanglement of 
a transformed separable state \ip') = U\ip sep ). For in- 
stance, the average linear entropy of a random product 
state transformed by a bi-partite unitary gate leads to its 
entangling power, introduced by Zanardi et al. , and 
later investigated for various models in [Tol - fT6| . An al- 
ternative approach to the problem of non-locality based 
on the mimimal Frobcnius distance of an analyzed global 
unitary matrix to the closest local gate was recently dis- 
cussed in , where a relation to matrix product opera- 
tor formalism was established. 

A given unitary gate U is called a perfect entan- 
gler, if there exists a separable state transformed by U 
into a maximally entangled state [ll]. Another class of 
maximally entangling unitary gates was characterized in 
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Two unitary gates are called locally equivalent if they 
coincide up to local transformations. The general prob- 
lem of finding necessary and sufficient conditions for local 
equivalence remains open. The full answer is known in 
the simplest case of a two-qubit system, since a canonical 
form of a unitary gate of size four is established and any 
gate can be uniquely described by a three-components 
vector called information content YA 121K23I . 
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In this work we analyze properties of a 'typical' quan- 
tum gate. In other words, we are going to average quan- 
tities characterizing each gate over a unique, unitarily 
invariant, Haar measure on the space of unitary matri- 
ces. An ensemble generated according to this measure 
is often referred to as circular unitary ensemble (CUE), 
the spectrum belongs to the unit circle and the ensem- 
ble is invariant with respect to unitary transformation. 
An exemplary algorithm of generating unitary matrices 
from this ensemble was discussed in j24j . To study the 
set of quantum gates we found it useful to define a spe- 
cial circular unitary ensemble (SCUE) containing special 
matrices with determinant equal to unity. 

One of the main result of this paper consists in deriv- 
ing the probability distribution P(a) for the information 
content of a random two-qubit gate, induced by the Haar 
measure on the unitary group. As an application of this 
result we compute the relative volume of the set of perfect 
entanglers with respect to this natural measure. Basing 
on numerical results performed for unitary gates of larger 
dimensionalities we are in position to predict asymptotic 
behavior of the average entanglement entropy of a ran- 
dom unitary gate. 

Furthermore, we analyze the class of unistochastic op- 
erations, introduced in [25|, which can be described by 
a coupling with an A-/-dimcnsional environment initially 
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in the maximally mixed state, 

p' = $ u (p) = Tr env [U(p®^)tf} (1) 

The partial trace is performed over the environment de- 
scribed in the Hilbert space %. If the dimension N of 
the principal system and the dimension M of the ancil- 
lary system are equal, the map is called unistochastic, 
while it is called k-unistochastic if M = kN. Thus a 
unistochastic map tyjj is determined by a unitary matrix 
U of size ./V 2 , while any /c-unistochastic map is given by 
a matrix of size TV fe+1 . 

By construction any unistochastic map is bistochas- 
tic, since the maximally mixed state is preserved, 
$u(l/N) = Tr B (ljv2)/iV 2 =1 N /N. On the other hand, 
the converse is not true, and in this paper we deter- 
mine sufficient and necessary condition for a one-qubit 
bistochastic map to be unistochastic. 

The name of this class of maps is related to the clas- 
sical case, in which probability vectors are transformed 
by stochastic matrices. The matrix is called bistochastic 
(or doubly stochastic), if it preserves the uniform (max- 
imally mixed) probability vector. A bistochastic matrix 
B is called unistochastic if there exists an unitary V of 
the same size such that Bij = |T4,| 2 . For N = 2 all 
bistochastic matrices are unistochastic, but this is not 
the case already for N = 3. For higher dimensions the 
problem finding necessary and sufficient conditions for 
unistochasticity remains open [26j . 

In analogy to the classical case, any quantum map de- 
termined as in ([!} by an orthogonal matrix will be called 
orthostochastic. In [2|| it was shown that for any unis- 
tochastic map the spectrum of the corresponding dynam- 
ical matrix is given by the Schmidt coefficients of the 
unitary matrix U treated as an element of the composite 
Hilbcrt-Schmidt space. This implies that the entropy of 
such an operation S(^u) is equal to the entanglement 
entropy of the unitary matrix. In other words, a link be- 
tween non-local properties of a unitary gate acting on a 
bipartite system and the decoherence induced by an as- 
sociated unistochastic map acting on a single system can 
be established. 

It is appropriate to mention that the class of quantum 
maps for which the system is coupled with the environ- 
ment in the maximally mixed state was already investi- 
gated in the literature. Such maps were discussed in the 
context of quantum information processing (27l . |28| , and, 
under the name 'noisy maps', while studying reversible 
transformations from pure to mixed states [29|]. More- 
over, Haagerup and Musat [3(| analyzed properties of 
a related class of factorizable quantum maps introduced 
in |3lj . In fact fc-unistochastic operations coincide with 
a subset of these quantum maps called exactly factoriz- 
able. In general this set is not convex, and its convex hull 
defines a larger set of maps called strongly factorizable. 

Note that a given unitary matrix U of a composite 
dimension d = N 2 may play very different roles in the 
theory of quantum information. Let us specify here three 
most natural applications: 



a) U € U{N 2 ) describes a quantum gate acting on 
a N x N bi-partite system, and its operator Schmidt 
decomposition characterizes the non-local properties 0, 

b) U € U{N 2 ) determines by Eq. ([]} a unistochastic 
quantum operation, p' = <&jj(p), acting on an iV-level 
system [25j . 

c) U G U (N 2 ) defines a maximally entangled state of a 
composite, N 2 x N 2 system (HHH, as \tp) = (f7<8>I) | V+> 
where \ip+) = jf £^Li \ 

The paper is organized as follows. Unistochastic oper- 
ations are analyzed in section II. As any one-qubit unis- 
tochastic map is determined by a unitary matrix of order 
four, we analyze in section III the set of all two-qubit uni- 
tary gates. This allows us to characterize the set U2 of 
one-qubit unistochastic maps, which forms a non-convex 
subset of the tetrahedron of bistochastic maps spanned 
by three Pauli matrices and the identity map. In section 
IV the ensemble of random two-qubit quantum gates is 
described. It is based on circular unitary ensemble of 
special unitary matrices of size four. The probability 
distribution of purity (nonlocality) for this ensemble is 
computed. Furthermore, we derive the probability that 
a generic gate belongs to the class of perfect entanglers, 
so it can transform a product state into a maximally en- 
tangled Bell-like state. Unitary gates acting on TV X N 
systems are investigated in Section V. For completeness, 
some basic properties of the operator Schmidt decompo- 
sition and related algebra of reshuffling of a matrix are 
reviewed in Appendix A. 

II. UNISTOCHASTIC MAPS 

Any unitary matrix U of size N 2 describes a unitary 
gate acting on the bi-partite system. Alternatively it 
may be used to define a unistochastic map [25| acting on 
a single system of size TV according to eq. (pj . In other 
words, the principal system is coupled to the ancilla of the 
same size, M = N, prepared initially in the maximally 
mixed state. Unless the gate U is local so that U = 
U a ®Ub, the partial trace leads to a non- unitary evolution 
of the density matrix p. 

An y su ch discrete map can be written in the Kraus 
form [H 

k 

p' = $(p)=J2AiP4 ■ (2) 

i=l 

To preserve the trace, Trp' =Trp = 1, the Kraus opera- 
tors need to satisfy the completeness relation 

k 

E4^ = i- (3) 

A trace preserving map $ written in the form ([2]) is called 
stochastic. It can be represented by the dynamical ma- 
trix D$ = TV(<f> g) I)\ip + )(ip + \, also called Choi matrix. 
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To perform the partial trace over the environment in 
the definition ([T]) let us apply the operator Schmidt de- 
composition of U recalled in (|A3| . making use the nota- 
tion introduced in Appendix the map $ can be rewritten 

as 



p' = $up = Tr c 

N 2 N 2 



= 1 3=1 



(4) 



The standard Kraus form is obtained by rescaling the 
operators, Aj = y/ Ai/NB' i: where B' i7 which arise by re- 
shaping the eigenvectors of (U R )^U R according to (|A4|) . 
Reshaped Kraus operators are eigenvectors of the hermi- 
tian dynamical matrix D mi i of size N 2 which determines 

}35j the map 



$(/o) so 



3 mu — Dmn Pnv 



(5) 



Taking into account an appropriate normalization we ob- 
tain the dynamical matrix corresponding to the unis- 
tochastic map 



D* v = ±(U R )i U R 



(6) 



Hence the rescaled Schmidt coefficients Aj/AT of any uni- 
tary matrix U treated as an element of the composite 
space of matrices of size A^ 2 provide the spectrum of the 
dynamical matrix D representing the action of the unis- 
tochastic map *5>u- 

The matrix D is normalized according to TrZ? = 

2i=i = Therefore the entanglement entropy 
(|A8|) characterizing the non-local properties of U is equal 
to the entropy of an operation &u j25(. For instance, the 
entropy vanishes for any local U which induces an unitary 
$[/ and is equal to 2 mA' for the Fourier matrix (|A10[) , 
corresponding the the maximally depolarizing channel, 
*f(p) = 1/N. 

The environmental representation (fTJ may be gen- 
eralized by allowing a larger size of the environment. 
Their physical motivation is simple: not knowing any- 
thing about the environment (apart from its dimen- 
sionality), one assumes that it is initially in the maxi- 
mally mixed state. In particular, one can define gener- 
alized, K -unistochastic maps determined by a unitary 
matrix U(N K+1 ), in which the environment of size N K 
is initially in the state 1 n k/N k . By definition, a 1- 
unistochastic map is unistochastic. 

A map $ is called bistochastic if it preserves the 
trace and keeps the maximally mixed state invariant, 
$(1/N) = 1/N. The definition (TTJ implies that both 
these conditions are satisfied, so any unistochastic map 
is bistochastic. In the following sections of this work we 
demonstrate that the converse is not true. 



In the simplest case of one qubit maps any bistochastic 
map is called a Pauli channel, since it can be brought by 
means of unitary rotations into the form (3(| 



P -> P 



2J A, a 



iPOl 



with 



£ Ai = 1 . (7) 



Here <7i denote Pauli matrices while ctq = 1. The Pauli 
matrices satisfy o~j = —iexp(iTT<jj/2), so the extreme 
points of the set of the bistochastic maps represent rota- 
tions of the Bloch ball around the corresponding axis by 
the angle 7T. 

Describing density matrices through their Bloch vec- 
tors, p = |l + t • a, we can write any bistochastic map 
p' = <f>/9 in the form 



tf 



(8) 



Here t denotes a real matrix of size 3 which we bring to a 
diagonal form by orthogonal transformations 0\ and O2 ■ 
As we permit only unitary rotations of the qubit repre- 
sented by orthogonal matrices from SO(3), (reflections 
are not allowed), some elements of the diagonal matrix 
•q may be negative. The elements of the diagonal matrix 
r\ are called the damping vector if = (j]x, 772, %), because 
the transformation ([5]) takes the Bloch ball to an ellipsoid 
with three axis given by ff. 

Writing down the supcropcrator $ of a Pauli channel 
and reshuffling this matrix according to (|A6|) we obtain 
the dynamical matrix D. For any bistochastic maps the 
dynamical matrix D splits into two blocks and its eigen- 
values are 



dx, 



-A 



f 

2' 
f 



1,4 



l 2,3 



-A 



2,3 



Vz ± {Tlx 
- Vz ± (Vx 



-Vy)} 



(9) 



Due to Choi theorem the map $ is completely positive if 
the dynamical matrix D is positive definite. This is the 
case if all eigenvalues d% = Aj/2 are not negative, which 
is true for the damping vector ff satisfying the Fujiwara- 
Algoet conditions [37[ 



(1±%) 2 > (m±V2) 



(10) 



These four inequalities assure that the corresponding 
positive map is CP. They define a regular tetrahe- 
dron whose extreme points are ff = (1, 1, 1), (I, —I, —I), 
(—1,1,-1), (—1,-1,1). The first point represents the 
identity operation while the three others correspond to 
unitary rotations by one of three Pauli matrices. As 
shown in section fill Bl there arc bistochastic maps given 
by ff inside the tetrahedron, for which the representation 
(H|) does not exist, so they are not unistochastic. 



III. TWO QUBIT UNITARY GATES 

Treating unitary matrices as quantum gates we need 
not to care about an overall phase, since physical states 
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differing by such a phase are identified. Discussing gates 
acting on a bipartite system we may thus fix this phase 
and restrict our attention to matrices pertaining to the 
special unitary group SU(N 2 ). 

In this section we are going to analyze the simplest 
case of two-qubit unitary gates. Hence we set N = 2 
and study matrices of the group SU (4) . A choice of two 
subspaces distinguishes a subgroup SU (2)® SU(2). From 
the physical point of view it corresponds to selection of 
two distinct subsystems. 




(f . }.°) 




Canonical form 



Consider unitary matrices U and V of size N 2 



N 2 



which act in the composite Hilbcrt space T-Ln®T~Ln- Two 
such matrices are locally equivalent, written U ~i oc V, 
if there exist local operations, Wa ® Wb and Wq <8> Wr>, 
such that V = (Wa ® W B )U{W C ®W D ). 

It is known, 0, HH, HU , that any unitary matrix U of 
size 4 is locally equivalent to some matrix of the following 
canonical form, 



V = exp^i^ a k <Jk ® cruj ■= exp(iH int ^ (11) 
fc=i 

where stand for the Pauli matrices, and aj~, are real 
for k = 1,2,3. Here Hint denotes a Hermitian operator 
of order four, called the interaction Hamiltonian of the 
two-qubit gate. 

The two-qubit gate V is periodic in each parameter 
afe with period it/2. Thus the space of parameters char- 
acterizing the gate has a structure of a 3-torus T 3 . How- 
ever, the coefficients a*, are not determined unambigu- 
ously. Due to the symmetries of the problem out of a 
cube spanned by three components of a one can distin- 
guish a smaller subset, called Weyl chamber, such that 
the correspondence between its points and the local or- 
bits is one-to-one [2lJ. For instance, it is possible to 
bring locally the Hamiltonian Hi nt into a form in which 



j > ot\ > a 2 > a 3 > or (12) 
i > a 2 > a 3 > 0. (13) 



7T VI 

2^ ai> 4' 2 



7T 

— — ai 



These restrictions imply that the Weyl chamber has a 
structure of a tetrahedron, which forms a 1/24-th part of 
the cube - see Fig la. 

The three-component vector a, called interaction con- 
tent of the gate, characterize the purely non-local inter- 
action Hamiltonian Hi nt . By construction two unitary 
gates are locally equivalent if and only if they are char- 
acterized by the same interaction content. 

For completeness we present a direct algorithm of find- 
ing the information content a of any two-qubit unitary 
gate U based on 

dEl. It relays on the fact from the 
group theory, SU(2) x SU(2) ~ SO{4), which implies 



FIG. 1: a) Any orbit of locally equivalent two-qubit unitary 
gates intersects the Weyl chamber forming a tetrahedron in 
the cube of vectors a representing the interaction content, 
b) The set of perfect entanglers forms a polytope inside the 
Weyl chamber with corners at points P, Q, N, M, CNOT and 
DCNOT [3 . Note that point P represents VSWAP. 



that any local operation from SU (2) x SU(2) forms in the 
so-called magic basis an orthogonal matrix O G SO (4) 
[38| . Therefore any unitary U written in this basis can 
be brought to the canonical form (|11[) by an orthogonal 
rotation, 



V = T UO 



(14) 



The algorithm works for any U € U(A) and consist of five 
steps: 

i) Find U' — Ue~ lx / 4 with \ equal to the phase of the 
detC/, such that W G SU{4); 

ii) Write it down in the magic basis, i.e. find W = 

/ -i -i \ 



MU'M\ where M = -j= 



contains 



1 
i i 
\ 1-10/ 
(row-wise) four Bell states forming the magic basis 

{-#+>, |0+),-#->,|v->}; 

iii) Compute WW T and find its spectrum, which we will 



write in the form {e 



-2i<5i 



-2iS 2 



-2iS 3 



-2iS 4 



}; 



iv) Find vector S by dividing eigenphases of WW T by 
minus two. Pay attention to the total phase: By con- 
struction - Si should be equal to zero, so if this is not 
the case replace b max by S max - n or 8 min by 5 min + tt, 



which corresponds to another choice of the signs in v e . 
v) Change variables to obtain the information content a, 

at = (61 + 82 - S 3 - 6i)JA = (81 + 6 2 )/2, 
a 2 = {S 1 -S 2 + S 3 -S 4 )/A={S 1 +S 3 )/2, (15) 
a 3 = (-8 1 + 8 2 + 8 3 -8 i )/4=(8 2 + 8 3 )/2. 



Alternatively, the vector 8 can be defined as the spec- 
trum of the Hamiltonian Hi nt entering the canonical form 
(fTTj) . The Hamiltonian H int is traceless and diagonal in 
the magic basis. Its four eigenvalues Si depend on the 
information content a in a linear way. An inverse of |1E 
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gives 

51 = oti + a 2 - as, 

5 2 = a\ - u 2 + a 3 , 

5 3 = -ui + «2 + ce 3 , 

54 = —ax — a.2 — a 3 , 

so one may easily switch between both representations. 

Note that the existence of the canonical form (fTTj) for 
N = 2 is due to the fact that the group SU(2) is ho- 
momorphic to 5*0(3). However, for higher dimensions 
SU(N) is homomorphic to a measure zero, (N 2 — 1) 
dimensional proper subset of the (iV 2 — 1)(7V 2 — 2)/2 
dimensional group SO(N 2 — 1), and no direct analogue 
of such a canonical form is known. 



if the absolute value of arguments of arc cosine are 
smaller then unity. This leads to the following restric- 
tions for the set U of damping vectors corresponding to 
unistochastic maps 

mm < m 

mm < m (21) 
mm < m- 

A bistochastic one-qubit map $ r y is unistochastic, if the 
vector fj satisfies the set of three conditions (f2~Tj) . so rela- 
tion (|20[) gives the interaction content a and the explicit 
form of the unitary matrix defining the map. 



B. Schmidt coefficients 



By definition (| A3|) the Schmidt vector A = A/4 is in- 
variant with respect to local unitary operations. To find 
its a relation with the information content a of an ar- 
bitrary unitary matrix of size N = 4 it suffices to take 
U in its canonical form (fTTj) . and find singular values of 
the reshuffled matrix U R . Their squares appear in the 
spectral decomposition of a positive, Hermitian matrix 

U can (U^y = AiI<g)I + A 2 cri®CTi+A3CT2®0-2+A4Cr3«>cr 3 

(17) 

Simple algebra gives the connection between local in- 
variants: 



(-1,-1,1) 
[0,0,0, 1] 




(1,1,1) 

[1, o, 0, 0] 



(a) 



V2 



(-1,1,-1) 
JO, 0, 1, 0] 



Ax = (1 + cos 2a 2 cos 2a 3 - 
A 2 = (1 + cos 2a 2 cos 2a 3 ■ 
A3 = (1 — cos 2a 2 cos 2a3 - 
A4 = (1 — cos 2a 2 cos 2«3 



cos 2ai cos 2ai3 - 

- cos 2ai cos 2a>3 

- cos 2ai cos 2«3 

- cos 2«i cos 2a3 



Introducing new variables 



cos 2ai cos 2a 2 ) 

- cos 2a± cos 2a 2 ) 

- cos 2«i cos 2a 2 ) 
h cos 2«i cos2a 2 ) 

(18) 



cos 2a 2 cos 2a3 
cos 2ai cos 2a3 
cos 2a\ cos 2a 2 



(19) 




LOCAL 
[1,0,0,0] 



we realize that relations (|18|) are equivalent to ((9]). Thus 
the coordinates fj used above represent the damping vec- 
tor which enters Eq. (|8]) and characterizes the corre- 
sponding unistochastic map Note that the vector A 
of Schmidt coefficients of U, determines the eigenvalues 
di = Ai/N of the dynamical matrix Djj associated with 
while Xi = Ai/N 2 describe the weights in the Pauli 
channel (JT)). 

To see which damping vectors fj correspond to unis- 
tochastic maps we have to describe the image of the Weyl 
chamber (fT3"j) with respect to transformation (fT9"|) . We 
may invert this relation 




2" m 

iarccos 2 ^ 22 - 

172 



(20) 



rarccos 



FIG. 2: a) The set U2 of one-qubit unistochastic maps forms 
a proper subset of the tetrahedron of bistochastic maps, b) 
The part of the Weyl chamber of non-equivalent unistochastic 
maps. 



The set IA 2 of one-qubit unistochastic maps forms a 
proper subset of the tetrahedron of bistochastic maps 
bounded by three parabolic hyperboloids 771 = 772^3, 
m = Vim j an d m = mm- Note that the set U 2 is not 
convex: it contains all four corners and six edges of the 
tetrahedron, but its center, the point fj = 0, belongs to 
the boundary of IA 2 - sec Fig 2. Thus we may model 
U 2 by pressing symmetrically four faces of a tetrahedron 
with rigid edges in such a way that they all touch in the 
center - sec Figs. 2 and 3. 

Let us now try to invert relations (TT9"|) and (|2"U|) to 
obtain the vector a as function of the components of the 
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FIG. 3: A model of the set U2 of one-qubit unistochastic 
maps which forms a non-convex subset of the tetrahedron 
of bistochastic maps. The model relays on the fact that the 
hyperboloids (|21[) can be ruled by straight lines. 



Schmidt vector, 



a 3 



jarccosv^i 



arccos 



(C 2 
W3 



where 
where 
where 



w 2 = 
w 3 = 



_ (Ai+A 4 -2)(Ai+A 3 -2) 
2(Ai+A 2 -2) 
(A 1 +A 4 -2)(A 1 +A 2 -2) 

2(Ai+A 3 -2) 
(Ai+A 2 -2)(Ai+A 3 -2) 



2(Ai+A 4 -2) 



(22) 

The choice of the sign of the square root determines 
the sign of a 3 . Constraints for unistochasticity (|21[) im- 
ply that \wi\ < 1 for i = 1,2,3 so the function arccos 
is well-defined. Alternatively, these inequalities provide 
constraints A has to satisfy, to represent a Schmidt vector 
of a N = 4 unitary matrix. The choice (fl"3"|) of the domain 
containing a implies that the elements of the Schmidt 
vector are ordered non-increasingly, Ai > A2 > A3 > A4. 

A local unitary operation (a = 0) is of rank one, while 
a generic two-qubit unitar y ga te has Schmidt rank 4. As 
noted by Diir and Cirac |22l | and Nielsen and al. 



there are no unitary matrices of Schmidt rank 3. Indeed, 
by setting A 4 to zero, equations (f!?2"j) impose A 3 = 
as well. Interestingly, for two-qutrit composite systems 
there exist unitary gates of size 9 with Schmidt rank equal 
to r for each r = 1, . . . , 9 [39l ]. 

Observe that expressions (|2"2"|) are ill defined if all com- 
ponents are equal, A* = (1,1, 1,1). Hence one may ex- 
pect that this degenerated Schmidt vector of the maximal 
entropy, S = 2 In 2, corresponds to different, non-locally 
equivalent gates. This fact, discussed in [221 ] . may be 
demonstrated with explicit examples. Let us define the 
following two-qubit unitary gates 



C^cnot 



10 
10 
1 
10 



C^SWAf 



f^DCNOT 



10 

10 
10 
1 



10 
1 
10 
10 



(23) 



Their names are related to applications in the theory 
of quantum information [l[. The gate CNOT performs 
the control not operation: Ucnot\o-, b) = \a,a®b) where 
a is the control bit, b is the target bit and the sum is 
understood modulo two. Defining a symmetric CNOT 
operation with the role of bits reversed, UcNOT>\a,b) = 
\a © b, a), one defines the double CNOT gate by the 
composition C/dcnot = C^cnotC^cnot'- The SWAP 
gate may be written as a product of three CNOT gates, 
%wap = C/cnot'C^cnotC^cnot', while its action on any 
two qubits, UswAp\&,b) = \b,a) explains its name. Char- 
acterization of non-local properties of these exemplary 
two-qubit gates is provided in Table 1, complementary 
to the data listed in [4(| • F° r comparison we include also 
the so-called i?-gate which interpolates between CNOT 
and DCNOT and is optimal to simulate an arbitrary two- 
qubit gate |4l| . 



The data presented in the table may be easily obtained 
for an arbitrary unitary matrix U: the spectrum of UU R 
gives the Schmidt vector A, the phases UU T (normalized 
in such a way that their sum vanishes) provide vector 
S, while the information content a follows from relations 

Note that the matrices Z7dcnot and [/swap are invari- 
ant with respect to reshuffling, so the reshuffled matrix is 
unitary and all four its singular values are equal to unity. 
Reshuffled Fourier matrix remains unitary, so it has the 
same singular values. Hence the Fourier matrix, the gates 
DCNOT and SWAP are characterized by the same, max- 
imally mixed Schmidt vector A* = (1,1,1,1), but they 



carry different information content, and thus are not lo- 
cally equivalent [22| • Making use of (|22l) we infer that any 
gate with the information content a = (7r/4, 7t/4, x) with 
an arbitrary x is characterized by the same maximally 
mixed Schmidt vector A*. 

For comparison we have provided the data for ma- 
trices representing y/ CNOT and VSWAP, which may 
be obtained by replacing the fragments of (|2"3"]l con- 

1 



taining the NOT gate U] 



not — 



1 



by U 



Vnot 



In these cases taking the square root of 



1 + i i - 1 
1 - i 1+i 

a gate corresponds to dividing its interaction content by 
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TABLE I: Nonlocal properties of two-qubit unitary quantum gates; s± — 2 ± y/2, t = 1 / %/2; perfect entanglers: No (N), Yes - 
inside the set (Y), Yes - at the boundary (B). 



Gates 


information 
content a 


Hamiltonian 
eigenvalues 5 


Schmidt 
vector A 


Schmidt damping 
rank vector ff 


perfect 
entangler 


local gate 


(0,0,0) 


(0,0,0,0) 


(4,0,0,0) 


1 


(1,1,1) 


N 


VCNOT 


1(1,0,0) 


1(1,1,-1,-1) 


( S+ ,s_,0,0) 


2 


(l,t,t) 


N 


CNOT 


1(2,0,0) 


f(2,2,-2,-2) 


(2,2,0,0) 


2 


(1,0,0) 


B 


B-gate 


1(2,1,0) 


f(3,l,-l,-3) 


|(3,3,1,1) 


4 


1(1,0,0) 


Y 


DCNOT 


f(2,2,0) 


1(4,0,0,-4) 


(1,1,1,1) 


4 


(0,0,0) 


B 


VSWAP 


1(1,1,1) 


1(1,1,1,-3) 


|(5,1,1,1) 


4 


1(1,1,1) 


B 


SWAP 


1(2,2,2) 


f(2,2,2,-6) 


(1,1,1,1) 


4 


(0,0,0) 


N 


Fourier 


f(2,2,-l) 


f(5,-l,-l,-3) 


(1,1,1,1) 


4 


(0,0,0) 


N 



two. a(U V swAp ) = ^a(U S wAp)- More generally, among 
several possibilities of taking a fc-th root of an unitary 
matrix U described by cv one can select an unitary matrix 
U 1 ^ such that its interaction content equals a/k. 

Let us mention here that the entangling power of in- 
teger roots of the SWAP gate were investigated in [42l ]. 
Another possibility to characterize nonlocal properties of 
unitary gates by so-called Frobcnius fidelity was recently 
investigated in [l7j ■ 



IV. SCOE(4) - ENSEMBLE OF TWO-QUBIT 
UNITARY GATES 



In this section we study random two-qubit unitary 
gates described by unitary matrices of size N = 4. An 
ensemble of unitary matrices generated according to the 
Haar measure is called [43j circular unitary ensemble 
(CUE). 

We aim to derive the probability distributions of local 
invariants of unitary matrices of size 4. Since quantum 
states are defined up to a global phase, we may restrict 
our attention to the set of special unitary matrices SU (4) 
with dct[7 = 1. An ensemble containing special unitary 
matrices with detU = 1 induced by the Haar measure 
will be called SCUE. 

To find the distribution P(S) we use magic basis and 
the representation given by Eq. (|14p . The Haar measure 
on CUE (4) induces in set of symmetric unitary matrices 
Y = UU T the measure 



d H Y = d H T YO . 



(24) 



invariant with respect to orthogonal similarity Y — > 
O t YO. A set endowed with such invariant volume is 
called circular orthogonal ensemble (COE), so the en- 
semble of symmetric unitary matrices with determinant 
equal to unity will be called SCOE. 

It is known [43j | that random symmetric unitary ma- 
trices of COE are characterized by the following joint 



distribution of eigenphases 6^ 

P(6i, 62,63,94) = MS (9! + e 2 + e 3 + e 4 ) 

IJ | e * e --e ie "| , 

l<m<n<4 

where M is the normalization constant. 
Integrating over 94 we obtain 



(25) 



P(6 l5 62,63) 
,i(-e 1 -e 2 -e 3 ) _ 



= TV' 



,i(-ei-e 2 -e 3 ) 



i(-ei-e 2 -e 3 ) 



,'C-) 2 



je 3 



e ,e 3 
„»©i l 



(26) 



The angles ir/2 > 6; > 7r/2, i :— 1, 2,3,4 are the eigen- 
values of Y = UU T , hence 6* = 28, [|. Thus 



Oil 

fl'i 
0.3 



1 



7 (Si + 5 2 



h - h 
h + <^3 - h) : 
S3 - Si 



5i) = - (9i + 9 2 ) (27) 



4 (9i + 



i.e. 



61 = 2(ai + a 2 — as) 
O2 = 2(ai — cc2 + Q3) 
93 = 2(— ai + Ci2 + 03). 



9 3 ) (28) 
+ 6 3 ), (29) 



(30) 

(31) 
(32) 



In this way we obtain one of the main results of this pa- 
per: the joint probability distribution for the interaction 
content vector cv, induced by the Haar measure on SU(A) 

P(ai,a 2 ,a 3 ) = N" |sin(2(ai +a 2 ))| |sin(2 (ai +a 3 ))| 
x |sin(2(a 2 +"3))! |sin(2(ai -a 2 ))| 
x |sin(2(a 1 - a 3 ))| |sin(2 (a 2 - a 3 ))| .(33) 

The resulting normalization constant M" can be easily 
calculated to be 2/7T by integrating over the cube — tt < 
®i < tt, i.e. -tt/2 < c 4 < tt/2, i = 1,2,3. 
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Making use of (|T9)) one can change variables and get 
the probability distribution of the damping vector ff, 



V1V2 7j3?ji 




V3V1 V2V3 




V2V3 V1V2 


V3 m 




m m 







1 



Vl 



1 



V3 



1 



Vl 1 12 

V3 



(34) 

defined for vectors ff satisfying constraints (f^Tj) . This 
expression shows that the density is concentrated in 
vicinity of the boundary of the set U 2 inscribed inside 
the tetrahedron of one-qubit bistochastic maps. Alterna- 
tively, one may easily get the analytical expression P(X) 
for the density inside the Schmidt simplex induced by the 
Haar measure on SU (4) , but we found it more convenient 
to work in the ff representation and to use (|34p . 

This distribution may be applied to compute average 
values of various measures of nonlocality of random two- 
qubit unitary gate. Consider, for instance the purity r 
(locality) related to the linear entropy of the Schmidt 
vector 

r := \\ + Xl + A3 + \\ = \ (1 + vl + vl + vl) . (35) 

The average purity of a unitary matrix distributed ac- 
cording to the Haar measure on U(N 2 ) was computed by 
Zanardi Q 



(r)N = 



N 2 + l 7 



(36) 



so in the case of two-qubit gates this average reads (r) 2 = 
2/5. 

Making use of the joint distribution (fM]) we can go 
a step further and get the entire probability distribu- 
tion P(r) as a triple integral. Two integrals are easy 
to perform analytically, but the last integral had to be 
computed numerically. Results shown as a line in Fig. IIVI 
are compared with the Monte Carlo calculations in which 
10 5 random unitary matrices of size 4 where generated 
and the distribution of their locality collected into the 
histogram denoted by black dots. Results obtained with 
both methods agree well and show that a typical random 
unitary gate has low purity, so it is strongly non-local. 
There exits a large class of quantum gates with the lin- 
ear entropy r ~ 3/8, close to the value for which the 
analyzed probability distribution P(r) achieves its max- 
imal value. All these gates represent a 'generic' behavior 
in the entire set of two-qubit unitary quantum gates, but 
it is hardly possible to distinguish out of them a single 
gate with some special properties. 



V. THE SUBSET OF PERFECT ENTANGLERS 
OF TWO-QUBIT UNITARY GATES 

Among all two-qubit gates represented by unitary ma- 
trices of size 4, one may distinguish so-called perfect en- 
tanglcrs, which can produce from a product state a maxi- 
mally entangled, Bell-like state. A gate U belongs to this 




FIG. 4: Probability distribution of purity r of random uni- 
tary matrices distributed according to the Haar measure on 
5(7(4). Solid line represents numerical integration of distri- 
bution (|34[) . while black dots show the histogram obtained by 
a Monte-Carlo approach. 



class if its numerical range (for a normal matrix equal to 
the convex hull of the spectrum) includes the eigenvalue 
z = 18]. In terms of the parameters ot\, a 2 ,a 3 the 
condition reads (3 



7T 


7T 


7T 


7T 


4 ~ 




> 

- 8 





> a 3 > 0, 



(37) 



In the latter work the authors have shown that the 
relative volume of the set of perfect entanglers is equal 
to 1/2. However, they used the uniform measure in the 
3-dimensional space of the vectors representing the infor- 
mation content (or rather its relevant part, called Weyl 
chamber), which does not corresponds the natural Haar 
measure on the set of unitary matrices. 

In this section we redo the calculations using the mea- 
sure P(a±, a 2 , ct 3 ) induced by the Haar measure on U(N) 
and given by Eq. (|3"3"|) . 

As explained in Section Mil the set of noncquiva- 
lent gates is parametrized by 0.1,0.2,0.3 restricted by 
inequalities (p~3|) . i.e. belonging to the tetrahedron To 
with vertices O = (0,0,0), A = (tt/2,0,0), C/dcnot = 
(7r/4, 7r/4,0), and [/swap = (tt/4, 7t/4, tt/4). Hence its 
volume with respect to the measure (|33[) equals: 



tt/4 q 1 q 2 

V w = j dai J da 2 J da 3 P (ai, a 2 , a 3 ) 
000 

7r/2 7r/2— qi q 2 

+ J dai J da 2 J da 3 P (ati, a 2 , a 3 ) = gj. (38) 

tt/4 

The set of the perfect entanglers (PE) given by (f3"T)) 
is the convex polyhedron with vertices: L = (7r/4, 0,0), 
M = (37r/8,7r/8,0), [/dcnot, Q = (tt/8, tt/8, 0), 
N = (3tt/8, tt/8, tt/8), and P = (tt/8, tt/8, tt/8), and a 
straightforward calculation of its volume with respect to 
the measure (|33|) gives: 



(39) 
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Hence the relative volume of the perfect entanglers reads 
^ = f « 0.85 , (40) 

and gives the probability with which a typical two-qubit 
gate is a perfect entangler. This value, much larger than 
one half, is consistent with the earlier observations, that 
a generic unitary gate is highly non-local. 




FIG. 5: A sketch of the set of two-qubit unitary quantum 
gates: the set of perfect entanglers (PE) occupies approxi- 
mately 85% of the entire volume and contains the set SPE of 
measure zero. The complementary set of non perfect entan- 
glers (NPE) consists of two (connected) parts which include 
local gates (Loc) and SWAP operation, respectively. 

In the set of two-qubit gates one distinguishes also 
the set of special perfect entanglers (SPE), which can 
maximally entangle a set of four product states which 
form an orthogonal basis [44[ • Special perfect entanglers 
form a one dimensional set of gates interpolating between 
CNOT and DCNOT, so the natural measure of this set 
is equal to zero. The interpolating gate B, introduced in 
[4l| , is located inside of the set of PE and belongs as well 
to SPE — see Fig. [5j which shows a sketch of the set of 
two-qubit unitary gates. 



VI. N x N SYSTEMS & TWO QUNITS 
UNITARY GATES 

Studying non-local properties of a unitary gate U of 
an arbitrary size N x N oi higher dimensions we have to 
relay on its Schmidt decomposition, since for A > 3 no 
direct analogue of the canonical form (|II[) exists. The 
Schmidt vector, equal to the squared singular values of 
the reshuffled matrix U R , in several cases may be found 
analytically. 



A. Schmidt vectors for exemplary gates 

Analyzing, for instance, the general case of the swap 
operator [/swap acting now in the space of two quNits, 
SWAP | a, b) = \b,a), we see that U R WAP = E/swap- Hence 
the reshuffled matrix is unitary, so all A~ of its singular 
values are equal to unity and the entanglement entropy 
is maximal, <S( E/swap) = 2 In AT. Selecting one element 
out of N in each block of the first A rows of U it is 
straightforward to write down the other AM — 1 permu- 
tation matrices, which are invariant with respect to the 
reshuffling transformation, so their entanglement entropy 
is maximal. 

The same property is also characteristic of the Fourier 
matrix of order N 2 with entries F^i = jj exp(27rfcZ/A f2 ). 
Also in this case it is enough to see that the reshuffled 
matrix F R remains unitary, which implies S(F) = 2 In N . 
This result was earlier established in @,[39j]. The fact that 
the entropy of the SWAP gate and of the Fourier matrix 
are maximal, does not imply that both gates are locally 
equivalent. Up till now the general question which uni- 
tary gates of size KN > 4 are locally equivalent remains 
open. 

Among other two-quNits gates let us mention two pos- 
sible generalizations of the CNOT (XOR) gate 

U+\i,j) = \i,i®j) and U-\i,j) = \i,iej), (41) 

where i, j = 0, 1, • • • , N— 1 and both operations are taken 
modulo N. 

Since addition and subtraction of bits modulo two are 
equivalent, both gates coincide for N = 2 with the stan- 
dard CNOT gate defined in ((231). The gate £/_ may be 
called a 'controlled rotation', because the target bit j 
gets rotated depending on the value of the control bit i 
and (U-) d = 1. Inasmuch as the gate U+ is symmetric 
and forms an involution, (E/+) 2 = 1, it was called the 
generalized XOR gate by Alber et al. (4(|, who demon- 
strated that U + is capable of performing various tasks of 
quantum information processing and proposed a physi- 
cal realization of this gate based on non-linear optical 
elements. 

Both gates are permutation matrices with the block 
diagonal structure. In both cases all blocks of size A 
reshaped into vectors of length A~ become mutually or- 
thogonal, so the matrix U±(U±)^ is diagonal and consist 
of exactly N non-zero elements, each of them equal to N. 
This fact implies that the entanglement entropy of both 
gates is equal and reads 

S(U+) = S(U-) = In A . (42) 

B. Random gates and average values 

The Haar measure on the set of matrices of a compos- 
ite size MN determines a natural measure on the set of 
unitary quantum gates. A gate taken randomly with re- 
spect to this measure will be non-local with probability 
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equal to one, since the set of local gates of the tensor 
product structure, Uj_(N) <S> U 2 (M), forms only a zero 
measure in U(NM). On the other hand, one may ask 
to what extend a generic unitary gate is non-local. For- 
mulating this question more precisely, we will investigate 
mean values of entropies of entanglement averaged with 
respect to the Haar measure. 

The average purity (|36p implies that the average linear 

entropy E = 1-E£ A f rcads ( E )n = {N 2 -l)/(N 2 + l) 
Q • Interestingly, this average is equal to the mean purity 
of squared components of a random complex vector of 
size TV 2 , distributed according to the natural, unitarily 
invariant measure on the space of pure states. However, 
we will demonstrate that both distributions are different. 
In particular, our numerical results show that the average 
entropy of entanglement (|A8|) . behaves like 

(S)u~2lnN-±, (43) 

while the mean entropy of a random vector of the size 
./V 2 reads 

N 2 

(S)f = *(7V 2 + l)-*(2) = V - w 21n7V + l- 7 . (44) 

k=2 

In this formula, derived in [47j, ^(x) denotes the 
digamma function while 7 w 0.5772 is the Euler con- 
stant. 



5.5 
5 

Co" 4.5 
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N 

FIG. 6: Mean Renyi entropies of entanglement averaged over 
unitary matrices acting on symmetric N x N systems: (•) 
represents the mean Shannon entropy (S), while (□), (*) and 
x denote average Renyi entropy (S q ) for q equal to 2, 4 and 
8 respectively. Dimension N is plotted in the log scale, and 
solid lines represent the asymptotic behavior. 

To characterize statistical distribution of the distri- 
bution of Schmidt vectors of random unitary gates we 
computed the average moments (Y). A|) and the aver- 
age Renyi entropy (|A9|) . Our results show that the av- 
erage entropy for the N x N unitary gates behaves as 
(S q ) ~ 2 In N — c q - see Fig. [5] Since the dimension N is 
plotted in the logarithmic scale, relations (|44p correspond 
to straight lines in the graph. The values of the constants 



c q agree with predictions c\ = 1/2, c 2 = In 2 w 0.69, 
C4 = (In 14)/3 « 0.88 derived in [48| for the average Renyi 
entropy of mixed quantum states of size N 2 obtained 
by partial trace of random pure states of an extended 
system. This fact shows that a random positive matrix 
U R (U R y has properties of a Wishart random matrix: for 
a large matrix size the unitarity of U hardly influences 
statistical property of the corresponding reshuffled ma- 
trix U R , which can be treated as a typical non-Hermitian 
matrix from the Ginibre ensemble [49j , as their statistical 
properties coincide asymptotically. 

A generic quantum pure state \ip) of a bipartite N x N 
system is typically highly entangled, so the corresponding 
reduced density matrix, p = Trjv|^>)(?/>|, is highly mixed 
and its von Neumann entropy is close to the maximal 
value, In N. This known fact 0, [5(| is directly related 
to properties of random unitary operation: a generic 
unitary gate U acting on N x N system is highly non- 
local and the corresponding unistochastic operation tyy 
is strongly depolarizing. 

Above observations hold also for a general case of k- 
unistochastic channels or quantum gates acting on N x 
M systems. The Haar measure on the unitary group 
U{NM) induces a certain measure on the simplex Ajv2 
containing all possible Schmidt vectors of size N 2 . The 
larger size M of the auxiliary subsystem, the larger is the 
average Shannon entropy of the Schmidt vector A', equal 
to the average entangling power of the quantum gate. 
In the limit M — > 00 the average entanglement entropy 
tends to the maximal value, 21nA^. 



VII. CONCLUDING REMARKS 

In this paper we have analyzed the ensemble of unitary 
quantum gates distributed uniformly with respect to the 
Haar measure on the unitary group. As the overall phase 
of the matrix does not influence its action on a quantum 
state, we restricted our attention to unitary matrices with 
determinant equal to unity and considered matrices from 
special circular unitary ensemble (SCUE). 

A generic unitary gate U, acting on a bipartite N x N 
system, was shown to be strongly non-local. For instance, 
its entanglement entropy S(U), (also called Schmidt 
strength 0]), behaves as 2 In N — 1/2, which is also char- 
acteristic of random density matrices of size N 2 , dis- 
tributed according to the flat measure [48| . 

Any bi-partitc quantum gate U, represented by a uni- 
tary matrix of size TV 2 , determines a unistochastic oper- 
ation in which the initial state of size is coupled 
by U with the environment of the same size, prepared 
in a maximally mixed state (25| . The entanglement en- 
tropy of a gate S(U), equal to the entropy S($u) of the 
corresponding unistochastic channel is given by the 
Shannon entropy of the eigenvalues of U R (U R y /N 2 . The 
reshuffled matrix U R is non-Hermitian and for a Haar 
random unitary U , its statistical properties are shown to 
coincide with predictions of the Ginibre ensemble. Thus 
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the spectral density of a normalized, Wishart-like ma- 
trix U R (U R ) ji /N 2 , is asymptotically described by the 
Marchcnko-Pastur distribution. 

Analyzing in particular the set of two-qubit unitary 
gates we derived an explicit formula for the joint prob- 
ability density (|33p for the interaction content vector a 
characterizing a quantum gate, and the distribution (|34[) 
for the damping vector ff. Having these results at hand 
we were in position to analyze the subset of bi-partite 
quantum gates consisting of perfect entanglers - unitary 
gates capable to transform a separable state into a max- 
imally entangled Bell-like state. Using the known condi- 
tions for a two qubit gate to be a perfect entangler |4Hl8j|. 
we found that according to the natural Haar measure on 
the unitary group the set of perfect entanglers occupies 
approximately 85% of the entire volume of the space of 
unitary matrices of order four. This observation gives a 
concrete argument supporting the claim that local gates 
are rather exceptional, while a generic unitary gate is 
strongly non local. 

Although any two locally equivalent unitary matrices 
possess the same set of the Schmidt coefficients and gen- 
erate the same unistochastic map, the reverse is not true. 
For instance, as shown in Table 1, the following two-qubit 
gates: SWAP, DCNOT and the Fourier matrix are char- 
acterized by the uniform vector A of the Schmidt coeffi- 
cients so their entanglement entropy is equal to 2 In 2. Al- 
though these unitary matrices are characterized by differ- 
ent information content a, so they are not locally equiv- 
alent, these gates generate the same unistochastic map: 
the maximally depolarizing channel. 

Any unistochastic operation, determined by a unitary 
matrix acting on an extended space, is by construction 
bistochastic. However, not every bistochastic operation 
is unistochastic and can be obtained by a partial trace 
over the environment of the same size, initially prepared 
in the maximally mixed state. In the case of one qubit 
maps, there arc no unistochastic maps of rank three, rep- 
resented by a point belonging to the face of the tetrahe- 
dron of bistochastic maps, spanned by three Pauli matri- 
ces and identity - see Fig. 3. This is consistent with the 
known fact that there is no two-qubit unitary gates of 
Schmidt rank three 0, [22[ . For instance the symmetric 
Pauli channel, 

1 3 

P -> P = o a ' l P ai ( 45 ) 

6 i=i 

at the center of the face of the tetrahedron of bistochas- 
tic maps is located as far from the set Ui of one-qubit 
unistochastic operations as possible. 
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N N202-090-239 of Polish Ministry of Science and Higher 
Education and by SFB/Transrcgio-12 program financed 
by Deutsche Sondcrforschungs Gcmeinschaft is gratefully 
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Appendix A: Schmidt decomposition of a unitary 
operator 

1. Matrix algebra: reshaping 

Consider a rectangular matrix A^j, j = 1,...,N and 
i = 1, M. Equivalcntly, one may put its elements row 
after row into a vector 5k of size MA, 

Aij = 3(j_i)jv + j. (Al) 

For instance, a unitary matrix U of size M x M will be 
thus transformed into a vector u with M 2 components. 

Let "Hjy denote an A dimensional complex Hilbcrt 
space, and T-Lhs the corresponding N 2 dimensional 
Hilbert-Schmidt space of all linear operators acting on 
Hn- It is equipped with a scalar product A ■ B = 
(A\B) := trA'B, where A and B are arbitrary complex 
matrices of size N x A. 

Let V denote an auxiliary unitary matrix of size A 2 . 
Unitarity of V implies that its N 2 columns b m = Vi, n 
(or rows b m = V m i), m,i = l, ...N 2 reshaped into square 
A x A matrices B m as in (|A1|) . form an orthonormal 
basis in HhS, since (B m \B n ) := Tr B} n B n = S mn . Note 
that the matrices B m need not to be unitary. 

Consequently, the tensor products {B m ® B n ), m, n = 
1, . . . , A 2 , span an orthonormal basis in the composite 
Hilbert-Schmidt space Whs ® T~Lhs of size A 4 in which 
acts the original A 2 x A 2 unitary matrix V. 



2. Operator Schmidt decomposition 

Take a given unitary matrix U of size A 2 x A 2 we 
wish to investigate. It belongs to the composite Hilbert- 
Schmidt space Hhs ® Has and will be occasionally de- 
noted as \U). Let us write down its representation in the 
basis defined above, 

N 2 N 2 

\U) = Yl C ^n\B m ) ® \B n ), (A2) 

m— 1 n— 1 

where C mn = tr((_B m (X) B n )^U). The complex matrix C 
of size A 2 x A 2 needs not to be Hcrmitian nor normal. 
The Schmidt decomposition of \U) reads 

N 2 

\U) = J2V^\B' k )®\B' k '), (A3) 

k=l 

where ^/Kk are the singular values of C, (square roots of 
eigenvalues of the positive matrix C^C) and the basis is 
transformed by a local unitary transformation W a <S> Wb- 
Thus \B') =W a \B), and \B") = W b \B) where W a and 
Wb are matrices composed of eigenvectors of C^C and 
CC\ respectively. In a typical case of a non-degenerate 
spectrum of CC\ the Schmidt decomposition is unique 
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up to two unitary diagonal matrices, up to which the 
matrices of eigenvectors W a and Wb are determined. 

Note that the matrix C depends on the initial basis, 
{B m Cg) B n }^ in=1 , in which the analyzed matrix U is 
represented, while the Schmidt coefficients are ba- 
sis independent. Thus, for convenience we may analyze 
the special case in which the basis in T-Lrs is generated 
by the identity matrix, of size N 2 x N 2 . Then each of 
the N 2 basis matrices B n of size N x N has only one 
non vanishing element which equals unity. Let's denote 
B k = B m » = \m)(fi\, where k = N(m - 1) + fx. In this 
case the matrix of the coefficients C has a particularly 
simple form, C mfl = Tr(B m » ® B nv )U = U mn . 



3. Matrix algebra: reshuffling 

This particular reordering of a matrix, called reshuf- 
fling [231, wn l be denoted as U R := C. In general the 
notion of reshuffling is well defined if a matrix X acts on 
a composite Hilbert space, Hm <8> 'Hat- The symbol U R 
has a unique meaning if a concrete decomposition of the 
total dimension, L = MN, is specified. Similar reorder- 
ings of matrices were considered by Hill et al. [U [H[ 
while investigating CP maps and also in [53^57} to ana- 
lyze separability of mixed quantum states. 

The Schmidt coefficients of U are equal to the squared 
singular values of the reshuffled matrix, U R . Therefore 
the operator Schmidt decomposition (|A3[) of an arbitrary 
matrix X may be summarized by 



\B' k ) 
\B'l) 



singular values of X R 
reshaped eigenvectors of (X R yX 
reshaped eigenvectors of X R (X R y 

(A4) 

Note that the singular values of the reshuffled matrix, 
SV(X R ), are equal to square roots of eigenvalues of a pos- 
itive matrix (X R y X R . The initial basis is transformed 
by a local unitary transformation W a ® Wb, where W a and 
Wb arc matrices of eigenvectors of matrices {X R )^ X R 
and X R {X R y , respectively. If and only if the rank K of 
X R [X y is equal to one, the operator can be factorized 
into a product form, X = X\ <E> X2, where X\ = T^X 
and X 2 = TriX. 

To get a better feeling of the reshuffling transformation 
observe that reshaping each row of the initial matrix X 
of length N 2 according to (|A1[) into a submatrix of size 
N and placing it according to the lexicographical order 
block after block produces the reshuffled matrix X R . Let 
us illustrate this procedure for the simplest case N = 2, 
in which any row of the matrix X is reshaped into a 2 x 2 
matrix 



a 





'x u 


X12 


X21 


X22 


Y R — 


AT13 


Xu 


X23 


X24 


X31 


X32 


X41 


X42 




. A^ 33 


A^34 


X43 


X44 _ 



(A5) 



The operation of reshuffling could be defined in an alter- 
native way. Instead of reshaping the vectors of X into 
square matrices of size N one can reshape columns of 
X, which leads to another reshuffled matrix X R . In 
the four indices notation introduced above (Roman in- 
dices running from 1 to N correspond to the first subsys- 
tem, Greek indices to the second one), both operations 
of reshuffling take the form 



Xmn 



and 



Xmn '■— Xi/jj 
nv nn 



(A6) 



However, both reshuffled matrices are equivalent up to 
a certain permutation of rows and columns and trans- 
position, so the singular values of X R and X R are 
equal. It is easy to see that (X R ) R = X. In general, 
Af 3 elements of X do not change their position during 
the operation of reshuffling (these typeset boldface in 
(|A"5])); the other A^ 4 - A^ 3 elements do. The space of 
complex matrices with the reshuffling symmetry is thus 
2N 4 — 2(N 4 - N 3 ) = 2N 3 dimensional. Note that if X 
is Hermitian the reshuffled matrix X R needs not to be 
Hcrmitian. 



4. Entanglement entropy 

The Hilbert-Schmidt norm of any unitary matrix is 
\\U\\ = y/(U\U) = s/N. Computing the norm of the 
right hand side of (|A3j) we obtain 



k=l 



N z 



(A7) 



Thus the normalized vector A of the squared singular- 
values, Afc := Ak/N 2 , lives in the (A^ 2 — 1) dimensional 
simplex and may be interpreted as a probability vector. 
Iff there exists only one non-zero singular value, Ai = 1, 
then the unitary matrix has a product form, U = U a ®Ub- 
In such a case U is called a local gate and both operators 
obtained by partial tracing, U a = Tt(,{7 and Ut = Tr a C/ 
are unitary. 

In general, the vector of the Schmidt coefficients of an 
unitary matrix U acting on a composite N x N system 
conveys information concerning the non-local properties 
of U. To characterize quantitatively the distribution of 
A one uses the Shannon entropy, 



jv- 



S(U) :=S(X) = - ^A fc ln(A fe ) 



(A8) 



k=l 



called in this context entanglement entropy of U , (or 
Schmidt strength Q), and the generalized, Renyi en- 
tropics 



S q (U) := S q (X) 



In 



fe=i 



(A9) 
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which tend to S in the limit q 1. The entropy So, 
sometimes called Hartley entropy, is equal to InL, where 
L denotes the number of positive coefficients Xi, and is 
called Schmidt rank (or Schmidt number). 

The second order Renyi entropy S2 is closely related to 
the linear entropy E(U) = 1 — exp(— S2) used by Zanardi 

in Q. The quantity r = X)fcLi(^*:) 2 i s called purity while 
analyzing the vector of eigenvalues of an arbitrary density 
matrix p: the larger coefficient r, the more pure state. 
The maximal value, r = 1 is attained if and only if the 
state p is pure. In the present analysis of the unitary 
matrices, we shall stick to this name, although in this 
context r could be termed locality: the gate U is local 
if and only if 7' = 1. Another quantity called inverse 
participation ratio is usefull: R = l/r = exp(S , 2) varies 
from unity (local gates) to N 2 for the Fourier unitary 
matrices of size iV 2 defined by 

F^ 2) := ^exp^TrfcZ/iV 2 ) • (A10) 
To demonstrate this fact it is sufficient to notice that 



the reshuffled matrix F R remains unitary, so all its sin- 
gular values are equal to unity, hence the Schmidt vec- 
tor contains N 2 equal components and is maximally 
mixed. Some examples of unitary two-qubit gates and 
their Schmidt vector are collected in Table I. 



5. Local equivalence 

By virtue of the Schmidt decomposition if two gates are 
locally equivalent, their Schmidt coefficients (and thus 
the entanglement entropy) are equal. However, the op- 
posite is not true: there exist unitary gates with the same 
set of Schmidt coefficients, which are not locally equiva- 
lent p2j . Hence equality of Schmidt vectors characteriz- 
ing two unitary matrices of size N 2 is a necessary but not 
sufficient condition for their local equivalence. Sufficient 
conditions for local equivalence are known [H, Hlj only 
for N = 2. 
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